Prom random walk to single-file diffusion 
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We report an experimental study of diffusion in a quasi-one-dimensional (qfD) colloid suspension 
which behaves like a Tonks gas. The mean squared displacement as a function of time is described 
well with an ansatz encompassing a time regime that is both shorter and longer than the mean time 
between collisions. This ansatz asserts that the inverse mean squared displacement is the sum of the 
inverse mean squared displacement for short time normal diffusion (random walk) and the inverse 
mean squared displacement for asymptotic single-file diffusion (SFD) ( <x -}^ > = 2 ^ t + 2F *i/ 2 
where D is the qfD self-diffusion coefficient and F is the single-file fD mobility). The dependence 
of F on the concentration of the colloids agrees quantitatively with that derived for a hard rod 
model, which confirms for the first time the validity of the hard rod SFD theory. We also show that 
a recent SFD theory by Kollmann l] leads to the hard rod SFD theory for a Tonks gas. 

PACS numbers: 83.50.Ha, 82.70.Dd, 83.80.Hj 



The diffusion of particles in quasi-one-dimensional 
(qlD) pores and channels is a basic feature of ion trans- 
port in cell membranes, molecular motion in zeolites, and 
particle flows in microfluidic devices (see references in 
|2j). The unique feature that separates qlD diffusion 
from diffusion in higher dimensions is the geometric con- 
finement that forces the particles into a single file with a 
fixed spatial sequence. This confinement generates a self- 
diffusion mechanism that has different time dependences 
of the mean squared particle displacement in different 
time domains. 

For time intervals shorter than the time between par- 
ticle collisions, in the presence of a randomizing back- 
ground fluid (e.g. a colloid particle in a solvent), the 
probability density for the particle displacement is 
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where x is the displacement during time interval t = 
t\ — to, and D Q is the qlD self-diffusion coefficient. How- 
ever, the fixed spatial sequence of the particles severely 
restricts the possibility for large single particle displace- 
ments and, therefore, drastically reduces the diffusion 
rate at long time. An analytic description of ID dif- 
fusion in a system of hard rods with stochastic back- 
ground forces was first reported by Harris 0- Several 
other ID systems have been examined with a similar ap- 
proach 0, 0, IS 0, S IS [ZZ3| ; the results obtained converge 
to the same solution. For an infinite iD system the long 
time behavior of the probability density for displacement 
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where F is a iD mobility defined by 



F = F 



We denote the ID mobility of the hard rods by F HR . 
In Eq.© a is the particle length, p is the ID number 
density, I is the mean spacing between the particles, and 
t c = l 2 /2D is the mean time between collisions in the 
system. Equations and © draw a remarkably sim- 
ple picture of ID diffusion at long time: the self-diffusion 
process, determined by the width of the probability den- 
sity, is proportional to t 1 ! 2 (i.e. < x{t) 2 >~ i 1 / 2 ) , and 
the proportionality constant is determined by the short 
time individual particle dynamics. 

Recently Kollmann reported an analysis of the long 
time behavior of ID diffusion that is valid both for atomic 
and colloid systems [l|. For colloid systems he finds the 
asymptotic particle density function displayed in Eq.J2J 
with the ID mobility, denoted by F q , 
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where q, S(q), D c (q) are the momentum transfer, static 
structure factor, and the short time collective-diffusion 
coefficient in q-space, respectively. The small t, small q 
approximation for the dynamic structure factor, S(q,t) 
[Tl| , yields the relation 

S(q,t) = S(q)cxpl -?D e (q)t\ ■ (5) 

I ) t<t e ,<j<S47r/<T 

Kollmann's analysis predicts that the long time character 
of ID diffusion is determined by the short time collective 
dynamics of the system. 

Although theoretical analyses of ID diffusion have 
been reported for the past four decades, the first exper- 
imental studies were reported only in the past decade. 
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Studies of molecular diffusion in zeolites, and of col- 
loid particles confined in a channel lead to the result 
< x(t) 2 >^ i 1 / 2 at long time Very recently, 

Lutz, Kollmann and Bechinger [l4| reported the results of 
an experimental study of single- file diffusion in a strongly 
interacting colloid suspension. The ID mobility, deter- 
mined from < x(t) 2 >= 2Ft 1 / 2 at long time, agrees with 
that determined from Eq.(0} at short time, as predicted 
0. However, the ID mobility they find is only qualita- 
tively similar to F HR . 

The main difficulty encountered in the study of single- 
file diffusion is to obtain data at long time; this diffi- 
culty is most pronounced for low concentration samples. 
To obtain the required data one needs a long-lived ex- 
perimental system and stable instruments, such as those 
cleverly devised for the studies reported in references 
and Q. 

In this Letter, we report an experimental study of 
qlD diffusion in a weakly interacting colloid suspension 
confined in a narrow straight groove. We establish an 
ansatz that accurately approximates the qlD diffusion 
process from the short time region to the long time re- 
gion, thereby allowing us to study the long time single- 
file diffusion within a reasonable time frame (requiring a 
sample lifetime of ~1 hour), as well as diffusion in the 
cross-over time region. The experimentally determined 
qlD mobility of the system agrees quantitatively with 

F HR 

We note that Kollmann states that F q is not equivalent 
to F HR , and the experimental results in Q support this 
statement. However, we show that these two theories are 
equivalent when applied to a systeimsuch as ours, which 
obeys the Tonks equation of state [18| ■ 

Our experimental system consists of silica colloid 
spheres (density 2.2g/cm 3 ) suspended in water and con- 
fined in straight and narrow grooves. The grooves are 
printed on a polydimethysiloxane substrate from a mas- 
ter pattern fabricated lithographically on a Si wafer 
(Stanford Nanofabrication Facility). The small width of 
the groove (< 2a) prevents the spheres from passing one 
other, and gravity keeps them from escaping the groove. 
The spheres are very weakly attractive (< 0.4/cbT); the 
short-range attraction is derived from surface tension ef- 
fects 0. Digital video microscopy is used to directly 
track the time-dependent trajectories of the spheres along 
the groove (the motion transverse to the groove is very 
limited and, therefore, is not considered here). Details 
relevant to sample preparation and data analysis have 
been described elsewhere 0, . 

We have studied qlD diffusion at various colloid con- 
centrations, characterized by a line packing fraction r\ = 
pa = Na/L, where L is the length of the groove in the 
field of view, and N the number of spheres within L. We 
used two different silica colloid suspensions. For 7/=0.09, 
0.17, 0.20, 0.38, 0.57, and 0.70 the samples had silica 
spheres with diameter a\ — 1.58/xm±0.04/im in a groove 



that was 3/im±0.3/7m wide and deep, and 2mm long. For 
?7=0.73 and 0.986 we used silica spheres with diameter 
CT2 = 3.7/iTO ± 0.1/Ltm in a groove that was 5/im ± 0.1/im 
wide, 4/irn ± 0.5/im deep, and 10mm long. Care has 
been taken to assure that there were no blockages in the 
grooves. We used the large spheres for the higher con- 
centration samples because the small spheres could not 
be contained inside the grooves when r\ > 0.7. 
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FIG. 1: Typical probability density of particle displacement 
evolving with time (for 7/=0.57). The solid lines are fits of the 
data to Eq.©. 

The self-diffusion process is usually described by the 
self-part of the van Hove function, G s (x,t), which is the 
probability density for finding a particle at a point xq + x 
at time to + t given that it was at xq at to (G s (x,t) = 

(E£=i S i x + x ^o) ~ ^(to + t)])E3). Figure [J shows 
a typical G s (x,t) for our system, derived from time- 
dependent trajectories. The deviation of G s (x,t) from 

a Gaussian, characterized by a2(t) = (^ 3 <^(t] 2 > 2 — 1^ , is 
found to be negligible (0:2 (*) ^ 0.1). We therefore assume 
G s (x, t) to have the Gaussian form 

G s {x,t) = —= 1 cxpJ X ®, . I. (6) 

V ' ^2tt < xHt) > \ 2<x 2 (t)>j W 

Figure n a l so shows the Gaussian fits to G s (x,t). The 
mean squared displacement determined from the fitting 
is sensibly the same as that determined from < x(t) 2 >— 

i(EliMto + t)- Xl (t )} 2 ). 

Figure [3 shows < x(t) 2 > as a function of t at various 
concentrations, extracted from fitting G s (x,t) to Eq.©. 
Qualitatively, < x(t) 2 > is proportional to t at short 
time, changes smoothly to < x(t) 2 >^ f (7 <1) at later 
time, and reaches < x(t) 2 >^ t x l 2 at long time for the 
higher concentrations. Because of the expected trend 
in the behavior of < x(t) 2 > as a function of 7/, it is 
reasonable to postulate that with long enough time the 
low concentration samples will also exhibit < x(t) 2 >~ 
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t 1 / 2 . Accordingly, we use the following ansatz to describe 
< x(t) 2 > over the entire time range: 
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Equation JJJ leads to 
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By construction, Eq. JHJ satisfies both the short and 
long time limits, and it provides a characteristic cross- 
over time, t x = {F/D ) 2 . If F = F HR , then t x = t^ R = 
2t c /n (see Eq.©), so that for hard rods t x is, essentially, 
the mean time between collisions. The fits of < x(t) 2 > to 
Eq.lJSJ shown in Fig[3]indicate that Eq.JSJ) is a reasonable 
approximation for all time, and the fitting yields three 
pertinent parameters describing the qlD diffusion: the 
short time self-diffusion coefficient, D Q , the long time qlD 
mobility, F, and the cross-over time, t x . 
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FIG. 2: Mean squared displacement as a function of t at dif- 
ferent concentrations. Note that < x(t) 2 > for large spheres 
is scaled by the factor 02 The data (symbols) are shifted 
downward a factor of 3 from one another for clarity. The error 
bars are smaller than the symbols used. For t < Is the movies 
were grabbed at 30frames/s, and for t > Is the images were 
grabbed at Aframes/s and 5frames/s for small and large 
spheres, respectively (only a subset of the data are plotted 
for clarity). The solid lines are fits of the data to Eq.JSJ. 

When t] < 0.4 the fitted values for D a are D Q \ — 
0.11±0.005/zm 2 /s andL> o2 = 0.036±0.005/_tm 2 /s, respec- 
tively for the small and large spheres (the lower concen- 
tration data for large spheres are not shown here) . These 
values are, within the experimental precision, the same 
as those calculated for isolated colloids confined by the 
three walls of the groove and, therefore, are used to 
determine F HR in Eq. © and t x in Eq. © . At higher con- 
centrations the fitted self-diffusion coefficient is slightly 



smaller (~ 70% — 80%D o ), suggesting that hydrodynamic 
interaction between colloid particles comes into play even 
at the shortest time accessible in our experiment [13J, E2. 




FIG. 3: Quasi-ID mobility (solid circles) determined with 
the empirical expression (Eq.JSJ) as a function of concen- 
tration (F for the large colloids is scaled by the factor 
((Ji /<T2 ) \J D \ /D 2 ) . The solid line represents F HR . Other 



symbols represent F q determined from Eq.JlJ. 
zooms into the data at higher r\. 



The inset 



Figure |2| shows the fitted F as a function of 77. When 
0.17 < i] < 0.57, F — F HR within the experimen- 
tal precision. However, F « F HR /2 when rj — 0.09 
and F < 2F HR when r\ > 0.7. Using the fitted F 
and D we find t x = 205, 150, 107, 17, 4.5, 5.2, 3.8, 0.006s 
for r\ = 0.09,0.17,0.19,0.38,0.57,0.70,0.73,0.986, 
correspondingly to be compared with t RR — 
955,172,116,19,4.1,1.3,0.7,0.001s, respectively (note 
for large spheres t x and t RR are scaled by a factor 
{(Ji/cr 2 ) 2 (D o2 /D ol )). For rj > 0.7 we can force F = F HR 
by replacing the colloid diameter in F HR with a larger 
effective diameter. We speculate that at higher concen- 
tration the colloid-colloid interaction, though weak, must 
be accounted for. Since the first peak of the pair cor- 
relation function is at a separation slightly larger than 
the sphere diameter the effective sphere diameter is 
thereby increased. 

The accuracy of the fitted F and D Q depends on the 
range of t relative to t x R . If the range of t extends both 
to t < R and t > t% R we can extract F and D Q 
accurately from Eq.ljSJ); if not the values obtained are 
less accurate, as shown by the discrepancies between the 
fitted F and F HR for 77 = 0.09, and the fitted D a and 
expected D Q at higher rj. 

We now show that, as is to be expected, F HR = F q 
for a qlD system that is described by Tonks e quat ion of 
state /(l — pa) = pksT (f is the linear force) 18]. The 
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FIG. 4: The experimentally determined static structure fac- 
tor S(q) for small q as a function of concentration, compared 
with that derived from the equation of state for a Tonks gas. 



relative isothermal compressibility of a Tonks gas is 



Xt/xto = (1 - pvf = (1 - Vf = S(q) 



(9) 



q=0 



where \t — — xfi ^ or ^> and \To — {ksTp)^ 1 . Figure 
^Jshows S(q) for our system for small q (2n/q ^> a/2) as a 
function of r\; the agreement with Eq.© clearly indicates 
that our system behaves like a Tonks gas. The slight 
shift of experimental values of 5(g) to larger rj from that 
of a Tonks gas is consistent with the weak colloid-colloid 
attraction |19| . 

Kollmann's theory relates F q to the collective diffu- 
sion coefficient, D c (q), and the relative isothermal com- 
pressibility, S(0). Substituting 5(0) given in Eq.® and 
D c {q) = D H{q)/S{q) 11] (H(q) is the hydrodynamic 

factor) into Eq.gJ), we obtain F q = l^J^ = F HR , if 
H(q) — 1, i.e. if hydrodynamic interaction is negligible. 

We have calculated F q as follows. First, 
S(q, t) was determined from the trajectories us- 
ing S(q,t) 



^(p«(*)P-g(0))> where Pq{t) 

i rr /'J.M rm. CI 



^7= fdxexp(-iqx) J2k=l S l x ~ x k(t)]- Tnen S(q,t) 
was fitted to the short time approximation (Eq.JSJ)) to 
extract D c (q) at small q (q <Si Aw /a). Within the short 
time range 0.03s < t < Is, S(q,t) is well described by 
Eq.© except for the case r\ = 0.986. Finally, F q was 
calculated using Eq.|@} for all the concentrations except 
r] = 0.986. As shown in Figdl F q agrees with F HR 
within the experimental precision. The data in Fig|3| 
also show that F q depends on q, which we attribute 
to hydrodynamic interaction in the system. In a qlD 
system hydrodynamic interaction is screened on the 
length scale of the channel width, so it can be treated 
as generating a pair-interaction [17J. Then the effect of 
H{q) on F q is not significant. A full discussion of H{q) 



in the qlD system will be published separately. 

It is worth noting that < x(t) 2 > in Eq.lJSJ is the width 
of a Gaussian which is the product of the short time 
probability density (Eq.QJ) and the long time probability 
density (Eq.J2Jl). The success of the approximation given 
in Eq. © suggests that the van Hove function displayed in 
Eq. @ is valid for all time for our system. It is further im- 
plied that the randomizing background that determines 
the short time behavior and the correlated motion that 
determines the long time single-file diffusion are coexist- 
ing independent processes with time dependent weights. 
For t <S t x and t ^S> t x the system exhibits normal diffu- 
sion and single-file diffusion, respectively. However, for 
t ~ t x , the motion of a particle in ID is hindered by its 
neighbors and the short time displacement distribution 
is modified by the long time distribution. 
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